Dissipation-driven phase transition in two-dimensional Josephson arrays 
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We analyze the interplay of dissipative and quantum effects in the proximity of a quantum phase 
transition. The prototypical system is a resistively shunted two-dimensional Josephson junction 
array, studied by means of an advanced Fourier path-integral Monte Carlo algorithm. The reentrant 
superconducting-to-normal phase transition driven by quantum fluctuations, recently discovered in 
the limit of infinite shunt resistance, persists for moderate dissipation strength but disappears in 
the limit of small resistance. For large quantum coupling our numerical results show that, beyond 
a critical dissipation strength, the superconducting phase is always stabilized at sufficiently low 
temperature. Our phase diagram explains recent experimental findings. 



Dissipation due to the coupling with the surrounding 
environment is an unavoidable effect accompanying 
the operation of any microscopic or mesoscopic quantum 
device. The knowledge of the influence of dissipative ef- 
fects on quantum coherence and quantum phase transi- 
tions (QPT) is therefore essential to assess the reliability 
of such devices in performing tasks which strongly de- 
pend on the possibility to maintain entanglement (phase 
coherence), e.g., in quantum computation. Among the 
quantum devices that have already found wide applica- 
tion, many are based on a collection of regularly arranged 
or single small Josephson junctions HOJ. These are also 
among the candidates for the physical implementation of 
the so-called q-bits Q. 

Josephson junction arrays (JJA), are prototypical sys- 
tems displaying a quantum phase transition with a con- 
trol parameter tuning the strength of quantum fluctu- 
ations. This has become progressively clear after it 
was pointed out in the late 70 's that the charging en- 
ergy of Josephson-coupled superconducting grains could 
lead to the quenching of the collective superconducting 
phase. Since then, several studies [1 HI have been de- 
voted to characterizing the superconductor-normal (SN) 
transition in JJA. Among the systems studied, the two- 
dimensional (2D) ones are the most interesting as no true 
long-range order is possible at finite temperature while a 
genuine QPT occurs at T — 5j; moreover, 2D samples 
can be fabricated in a controlled way and experimentally 
characterized HHQ. 

The main effect of dissipation in JJA is that of quench- 
ing the quantum fluctuations of the phase variables 
(thereby enhancing those of the conjugate charges) thus 
stabilizing the S phase HE El El EI El El- This has 
had ingenious experimental confirmations, e.g., in a JJA 
coupled with a 2D electron gas substrate that allows one 



to tune the dissipation strength [jj , as well as in identical 
JJA with different built-in Cr shunt resistors Q . 

In this paper, we use a Fourier path-integral Monte 
Carlo (PIMC) approach to analyze the competition be- 
tween dissipation and quantum fluctuations in JJA, for 
strong quantum coupling, and in proximity of the QPT. 
In particular, we study the phase diagram as a function 
of temperature, quantum coupling and dissipation. We 
show how a large enough dissipation leads to the dis- 
appearance of the zero-tempcrature QPT, as well as of 
the reentrant low-temperature behavior displayed in the 
limit of large shunt resistance E3 ■ The phase diagram we 
obtain is in remarkable agreement with the experimen- 
tal one. Concerning the reentrance, we can explain the 
two distinct behaviors found in experiments; indeed, a 
non-monotonic (i.e., reentrant) low-temperature behav- 
ior of the array resistance has been observed in unshunted 
JJA [||, but not in shunted ones 

JJA are essentially described by the quantum XY 
model, whose coordinates and momenta correspond to 
the wave-function phases tfi and the net Cooper-pair 
number (charge) hi, respectively, of superconducting is- 
lands arranged on a lattice: 

ij id 

where [<pi,hj]=i5ij and d runs over the z nearest- 
neighbor displacements. The Josephson energy E 5 sets 
the energy scale, making it convenient to use the dimen- 
sionless temperature t — T/E^, while the charging en- 
ergy involves the capacitance matrix C'ij — Cly , with 

T^-* = (z Sij - J2d^iJ+d) + V $ij including the mutual 
capacitance C and the self capacitance Co =T]C ;in the 
experimental samples C is typically dominant i-e., 
77 <C 1. We consider here a square-lattice, z = 4. 
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FIG. 1: Phase diagram for the square-lattice JJA with 
77 = 0.01 for increasing values of the dissipation strength 
■j = Rq/Rs = (diamonds), 0.15, 0.25, 0.4, 0.45, and 0.5 (cir- 
cles, from left to right). Inset: Trotter- and size-extrapolated 
helicity modulus T(g,t) for 7 = 0.15, at £ = 0.2 and £ = 0.05. 



In the classical limit (large C) the charging term is 
thermodynamically irrelevant and the SN transition is 
a Berezinskii-Kosterlitz-Thouless (BKT) one ^tJj lead- 
ing from a high-temperature disordered phase to a low- 
temperature quasi-ordered phase at t„ =0.892 The 
superconducting phase is weakened |3j by the quantum 
fluctuations of the phase when the characteristic charg- 
ing energy E c — (2e) 2 /2C is comparable to the Joseph- 
son energy P,, i.e., when the quantum coupling constant 
g = y/E^jE~, is of the order of unity . The coupling 
parameter g can be varied in the fabrication of JJA, and 
its further increase can finally drive the zero-T system 
through a QPT at a critical value g c , numerically |2(ij 
estimated ~3.4. Near the QPT the system is charac- 
terized by a sharp enhancement of anharmonic quantum 
fluctuations (lif. 

In this paper we consider normal Ohmic shunt resistors 
R s as the source of dissipation, whose strength is mea- 
sured by the dimensionless parameter j = Rq/R s , with 
the resistance quantum RQ = h/(2e) 2 . It is well known 
that dissipation enhances the S phase, i.e., t c (g, r y) in- 
creases with 7, but reliable results could be obtained only 
in the low-coupling regime Mean field, renormal- 

ization group, and variational approaches 
predict the existence of a critical value ■y c —2/z=l/2 
above which the QPT disappears and for any value of g 
the system is in the S phase at sufficiently low temper- 
ature. Moreover, it is not clear whether the reentrant 
behavior (with the N phase reappearing at lower T) ob- 
served at 7 = in the proximity of g c Jllj disappears 
before or at the critical value 7 C . As these phenomena in- 
volve the strong coupling region where approximate the- 
ories give contradictory answers |2lj | , accurate numerical 
data are required for a real understanding, as well as for 
the interpretation of the experimental data. 

To this purpose we employ an efficient Fourier PIMC 



technique recently introduced 0, |2^ ■ This starts from 
the path-integral for the effective partition function Z = 
§ Vcpcxp { — S[<p] — [<£>]}, where the action S[<p] cor- 
responds to the Hamiltonian JJ) and S 1 [ip] is the bilocal 
action 0, 0, |2^| , namely 
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7 r ^ 1 

S[(f]=J du ^ ipi(u) (pj[u)-~22cos(p id (i 



^ 

i/t 



, (2) 



S^tp] = ^Jdudu' k(u-u') [(pid(u)-(pid(u')] 2 . (3) 



Here u <E [0, 1/t] is the (dimensionless) 'imaginary time', 
l Pid = l Pi~ < Pi+d: an d K(it)=7t 2 /8 sin 2 (7rfu) is the dissi- 
pative Ohmic kernel. Using Fourier-Matsubara variables 
it is easily shown that the influence action turns into 
a local form, suggesting that dissipation can be easier 
dealt with in Fourier space. Therefore, at variance with 
the standard PIMC algorithm, that samples the variables 
{tfit = <pi{t I Pi) , 1 = 1, P} after discretization of the 
interval u € [0, 1/t] in P slices of size 1/Pt (P being 
the Trotter number), we proposed to sample the P 
Fourier components of <pn. Choosing P — 2M+1 one can 
write 



M 



fit = Pi + 2 ^2 



k=l 
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Pik 
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P 
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■Pik 



sin— J . (4) 
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The 2M components { <p ik , ip ik } and the zero- frequency 
component (fi are sampled by the Metropolis algorithm. 
MC autocorrelation times can be significantly reduced by 
alternating Metropolis moves with microcanonical over- 
relaxed ones |lfil l24j. However, the real improvement in 
simulation efficiency arises from the fact that the move 
amplitudes can be independently chosen and dynamically 
adjusted for each Fourier component, thus correctly sam- 
pling also strongly fluctuating paths H^Ell- Eventually, 
the finite- P overall action S + S z reads 
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with the 'kinetic' matrix Tijk—^r- sin 2 ^ 
The last term in Eq. (0 containing the Josephson in- 
teraction is understood to be expressed using the expan- 
sion note that it does not burden the simulation as 
the {ifie} are stored and just updated in the component 
to be moved, at variance with the i nteg ral appearing in 
standard Fourier PIMC algorithms [26j. We performed 
an extensive set of simulations involving L x L square lat- 
tices with linear sizes up to L = 48 and Trotter numbers 
up to P = 201, setting 77 = 0.01. 

Our main results concern the phase diagram in the 
(i, g) plane for different values of the dissipation strength 
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FIG. 2: Temperature behavior of the helicity modulus T(£), 
on the 8x8 lattice and P —> oo, for g = 3.4 and different dissi- 
pation strengths, showing that rising 7 the reentrant behavior 
disappears and T(t) becomes monotonic for 7 > 0.2. 



7, and are summarized in Fig.H the data points were ob- 
tained for P— 101 as in Ref. 16 by fitting the finite-size 
scaling relation for the helicity modulus (or stiffness) per 
island T = (L 2 E,)- 1 [d 2 F/dq 2 ] q=0 , defined as the re- 
sponse of the free energy F(q) under twisting the bound- 
ary conditions as ipi — » ipi + qi. 

As shown in Fig. ^ in addition to generally stabilize 
the S phase, dissipation hinders the mechanism caus- 
ing the reentrance. However, the reentrance persists for 
small values of 7 and a finite value 7 > 0.2 is necessary 
to restore a monotonic critical line. The persistence of a 
reentrant normal phase for 7 = 0.15 is illustrated in the 
inset of Fig. ^ where the thermodynamic helicity mod- 
ulus, as obtained through a systematic finite-size scaling 
analysis, is plotted as a function of the quantum cou- 
pling: at g — 3.5 it appears that T keeps a finite value 
for f = 0.20, while it vanishes for £ = 0.05, thus signaling 
that the established phase coherence disappears again at 
low temperature. 

It is interesting to see how the temperature behavior 
of T(t) changes with 7 in the region of the reentrance. 
In Fig. 12 several finite-lattice data for g — 3.4 show that 
rising 7 removes the reentrance to disorder at low tem- 
perature, as if g were decreased \u> . From these data 
one can again roughly estimate that the threshold where 
the reentrance disappears is at 7 ~ 0.2. 

Further increasing 7 the critical line in Fig. ^ progres- 
sively changes its curvature and a flex appears for 7 > 0.4, 
signaling the incipient stabilization of a low-temperature 
S phase for any value of the quantum coupling, as we ar- 
gue below. In order to address this important issue, we 
study the dependence on 7 of the helicity modulus and 
related quantities at higher values of g, namely g = 5 and 
g = 10, and at a fixed low temperature t = 0.05. In Fig. 
we report the helicity modulus and the 'pure-quantum' 
spread of the phase difference between neighboring is- 
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FIG. 3: Finite-size data (L = 24) for T (top panel) and A* 
(bottom panel) vs 7, at t = 0.05 and g — 5, for different values 
of P. Inset: 1/P 2 extrapolation of A 2 . 
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FIG. 4: Finite-size data (L = 24) for T vs 7, at f = 0.05 and 
<7 = 10, for different values of P as in Fig. |H] Inset: T(t) 
for 7 = 0.55 (circles) and 7 = 0.50 (squares) for L = 12 and 
P = 201. The straight dotted line in the inset marks the 
universal- jump value 2t/n. 



and different Trotter numbers including the extrapola- 
tion to P— >oo. For increasing dissipation, T remains 
zero in an interval and then, around a crossover value 
7~0.3, it abruptly starts to increase; for slightly larger 
7 the critical curve t c (g, 7) is hence expected to cross the 
point (g — 5, t = 0.05). The crossover is clearly connected 
with the quenching of the pure-quantum phase fluctua- 
tions (lower panel of Fig. EJ- Just at the crossover, the 
PIMC data for A 2 , display the phenomenon of a markedly 
weaker convergence with P (see figure inset), signaling 
that fluctuations of high Matsubara modes become sig- 
nificant. We suggest that this reflects the proximity to a 
phase transition mainly driven by quantum fluctuations, 
which are in turn modulated by the dissipation. 

Since at g = 5 the system displays a BKT transition 
for 7^0.3, a much larger coupling is required to check 
whether 7 > 7 C = 1/2 in fact ensures low-i ordering. We 
therefore performed simulations for quantum coupling as 
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FIG. 5: Zero-i phase diagram. Squares: estimates of g c 
from the data reported in Fig.0 The dashed line extrapolates 
them to the expected behavior of g c (7). Full and open circles: 
N and S phase, respectively, as observed experimentally 0. 
The reentrant behavior occurs in the shadowed region. 



large as g — 10, as reported in Fig. 21 Again a crossover 
of T, signaling the proximity of the BKT critical line, 
shows up at 7~0.6: as t is still finite, this value is only 
an upper bound for 7 C . Indeed, as shown in the inset, 
the transition occurs at low t also for 7 = 0.5. Our re- 
sults are therefore consistent with j c = 1/2, in agreement 
with early predictions 0, Il3l fl4| . In addition, the reen- 
trance displayed at 7 = disappears well before 7 C = 1/2, 
so there is no evident connection between the two phe- 
nomena. 

The quantum phase transition occurs at the value 
g = g c {"f) where the critical temperature t c (g,j) van- 
ishes: assuming that 7 c = l/2 and using the data re- 
ported in Fig. ^ to estimate g c for some values of 7, we 
can draw the line of quantum critical points g c (7) that 
separates the S and N phase in the (7,3) plane. Fig. [S] 
evidences that the resulting zero-i phase diagram is in 
remarkable agreement with the experimental findings by 
Takahide et al. 

In conclusion, we have obtained the quantitative phase 
diagram of a resistively shunted 2D Josephson junction 
array, for several dissipation strengths 7. Our results in- 
dicate that the reentrant low-temperature normal phase, 
recently evidenced in the limit of infinite shunt resis- 
tance (7->0) H03, 

persists in a small but finite range 
of values of 7 < 0.2 (i.e., i? s >32M2) and 3.2< ff <3.5. 
This explains why the reentrance was not detected in 
the shunted J J A samples of Ref . 0, displaying resistances 
i? s <18kf2 and quantum coupling values (17 = 2.2, 4.2, 
5.5, 9.4) outside the above range. For 7 > 7 C = 1/2 we ob- 
serve the SN transition at finite t for very large coupling, 
so that our data validate the prediction that above 7 C 
the superconducting phase is always stabilized. Our ex- 
trapolated zero-temperature phase diagram explains the 
available experimental observations. 
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